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Abstract 

The main assumption of the model is that in soft processes mesons behave like 
systems made of valence quarks and an effective vacuum-like held. The 4-momentum 
of the latter represents the relativistic generalization of the potential energy. The 
electromagnetic form factors are expressed in terms of the overlap integral of the 
initial and hnal meson wave functions written under the form of Lorentz covariant 
distribution of quark momenta. The calculation is fully Lorentz covariant and the 
form factors of the charged mesons are normalized to unity at i=0. 
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In the lowest approximation of the standard model, the elastic electron-meson scat- 
tering is the result of the spontaneous one photon exchange between the electron and one 
of the elementary constituents of the meson. 

The Lorentz covariant, model independent parametrization of the hadronic matrix 
element entering the expression of the scattering amplitude is: 

< M(P')\ C/(+oo,0) J^(0) l/(0,-oo) |M(P) > 

= T" = k m f em (t) (P + Py (1) 

where km is the electric charge of the meson, U s (t,t') is the time translation operator 
describing the evolution of the meson under the action of strong forces, t = (P — P') 2 is 
the momentum transfer and the Lorentz invariant function f em {t) is the electromagnetic 
form factor which contains the whole information one can obtain on the meson structure 
from elastic electron scattering. 

Due to the local, elementary charcater of the electromagnetic current </^ m (0) the matrix 
element ([!]) is usually related to the probability of finding the recoiling quark in a meson 
with a momentum different from the initial one. It may be said that the form factor 
shows to what extent the initial system, where one of the quarks has been replaced by 
the recoiling quark, is a meson with another momentum. A natural consequence, strongly 
supported by the experimental data, is that the form factor decreases with t, because the 
larger is the momentum transfer, the harder is to incorporate the recoiling quark in a 
bound system. 

According to this picture the calculation of the form factors resorts to the evaluation 
of a kind of overlap integral and the main problem is to find a Lorentz covariant internal 
wave function for a system made out of independent constituents. This is not an easy 
matter, because as shown by the well known example of the Bethe Salpeter equation [l[ 
it is hard to solve this problem without introducing some unphysical degrees of freedom. 

Up to now, the most reliable results concerning the form factors -weak or electromagnetic- 
have been obtained by alternative methods, like, for instance, QCD sum rules 0, lattice 
calculations ||, chiral perturbation theory (CPT) [|J. There have also been proposed 
methods to calculate the overlap integral by making use of potential models [[J, 0, or of 
some particular reference frames, where the explicit form of the binding potential can be 
ignored §, §. 

The model we use in this paper is an effective model for hadrons as bound states of 
quarks, and, just like the chiral perturbation theory f9j, it is intended to complete the 
low energy picture of QCD. Its basic features do not follow from the symmetry properties 
of the underlying theory, as in the case of CPT, but from the general properties of the 
ground states. 

The fundamental assumption of the model is that at low energy hadrons reveal a stable 
strucure which looks as being made of valence quarks and of a vacuum-like field The 
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4-momentum carried by $ is the relativistic generalization of the potential energy in the 
quark system. 

The main reason for introducing this effective description is suggested by the examples 
taken from the relativistic field theory, where a bound state is not the instantaneous 
effect of an elementary process, but of an infinite series of elementary interactions 0. 
We conjecture therefore that the binding forces will never be "seen" on an instantaneous 
picture of a bound state because they are the result of a time average. Our specific 
assumption is that the effective component <3> represents the average over a time To of 
the elementary quantum fluctuations generatig the binding. The time To depends on the 
underlying dynamics and must be sufficiently long in order to assure a stable result. 

Another reason for introducing the effective field $ besides the valence quarks is that 
a system made only of on-mass-shell particles having a continuous distribution of relative 
momenta does not behave like a single particle because it does not have a definite mass 

§■ 

In agreement with these remarks, we work in momentum space where the mass shell 
constraints and the conservation laws can be easily expressed. 

The specific assumption of the model is that the generic form of a single meson state 
is PI 



\Mi(P)) = 

j^—f d 3 p q^ d*Q tfW(p + q + Q- PMp, q; Q) 

(2tt) 6 J e p e q 

xu(p)T M v(q) &(Q) a)(p)tf(q)\0) (2) 

where a + ,b + are the creation operators of the valence qq pair; u, v are Dirac spinors 
and is a Dirac matrix ensuring the relativistic coupling of the quark spins. The quark 
creation and annihilation operators satisfy canonical commutation relations and commute 
with $ + (Q), which represents the mean result of the elementary excitations responsible 
for the binding. Their total momentum is not subject to any mass shell constraint 
and, in some sense, it is just what one needs to be added to the quark momenta in order 
to obtain the real meson momentum. This is in agreement with our assumption that Q is 
the relativistic generalization of the potential energy. We shall suppose accordingly that 
Q is time like and, from stability reasons, Qo < 0. 

The internal function of the meson is the Lorentz invariant momentum distribution 
function ip(p, q; Q) which is supposed to be time independent, because it describes an 
equilibrium situation. This means that it does not change under the action of internal 
strong forces and hence the time evolution operators U s (t,t') in eq. ([D can be replaced 
by unity. The main role of ip is to ensure the single particle behaviour of the whole system, 
by cutting off the large relative momenta. 

In the evaluation of the matrix element (|I|) we shall use the cannonical commutation 
relations of the quark operators 
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m 

and the expression of the vacuum expectation value of the effective field which is 



{oi(fc),o}(g)} = {h(k),b}(q)} = (2tt) 3 ^ 5 t ^\k-q). (3) 



defined as follows 10 



(o| $(go $ + (g 2 ) |o> = 
(v^To)- 1 |rf 4 x e l ( Q2 - Qi )- XM = 

(2^) 4 (yT )- 1 ^(g 1 -Q 2 ) (4) 

where V is the volume of a large box and T is the characteristic time involved in the def- 
inition of the mean field $. It is important to remark that the definition (|]) is compatible 
with the norm of the vacuum state if one takes $(0) = 1. We notice also that the relation 
(||) has the character of a conservation law, just like the commutation relations (|3|), both 
of them being necessary for the fullfilement of the overall energy momentum conservation 
in the process. 

As a first test of the model we evaluate the norm of the single meson state (0) according 
to the usual procedure. The factor f T dX e % ( E ( p >~ E ( p coming from the 8^ functions 
in eqs. (|2|) and (f|) shall be put equal with T, because we assume that the incertitude in 
the meson mass is much smaller than T _1 . A short comment concerning this question 
will be given at the end. Observing that T is nothing else than the time involved in the 
definition of the effective field $ for a moving meson, we write it as T = J| T and get: 

( M(P') \M(P) ) = IE (2n) 3 5 i3) (P - P') J (5) 

where 

MV J e„ e„ 



X *W(p + q + Q- P)W(p, q; Q)\ 2 Tr [^±^ 75 ^-^ 75 j = I. 

This a remarkable result because it shows that the wave function of the many particle 
system representing the meson can be normalized like that of a single particle if the 
integral J converges. 

As a matter of consistency, we also remark the disappearance of the rather arbitrary 
time constant T from the expression ([5]) of the norm. 

We evaluate now the matrix element (0) proceeding in the same manner as before. By 
introducing the expression of the electromagnetic current written in terms of free quark 
fields 

J£n(x) = tto E nrfiW-fMn) (6) 
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between the meson states (|2[) and using the relations (||) and (^) to eliminate some integrals 
over the internal momenta, we obtain after a straightforward calculation: 



vt I d * Q lr It It 6i ^ p + q + Q ~ p) ^ q] Q) W + ^ 2)) (7) 

P Q fC 



A* 

2tt 



where 



= K 1 <J< 4 >(fc + g + Q-P') V>/(*,<7;Q) 



x Tr 



j 5 (k + mi)7 M (p + mi)7 5 (-g + m 2 ) 



t( 2 ) = « 2( 5( 4 )(p + & + Q-P')^/(P^;Q) 



x Tr 



7 5 (-A; + m2)7 M (-g + m2)75(j3 + mi) . (9) 



The two terms in (|7|) represent the contributions of the valence quarks, k is the mo- 
mentum of the quark after the absorbtion of the virtual photon, P' is the final meson 
momentum and ifij are the momentum distribution functions of the initial and final 
mesons respectively. 

In the next we shall work in the reference frame where the momenta of the initial 
and final mesons are P = (E, 0, 0, P) and P' = (E, 0, 0, — P) respectively and the 
electromagnetic form factor expresses as: 

femit) = — ^L—r %. (10) 

k m v4M 2 - t 

In this frame it is an easy matter to show that f em (0) = l. The demonstration makes 
use of 5^(k + q + Q) to eliminate the integrals over k in 7g and of the identity (p + 
m )lo{p + m) = 2e p (p + m) to reduce the number of projectors. Performing a similar 

(2) 

operation on T and proceeding like in the case of the norm, one gets 

T = 2M («! - k 2 ) J (11) 

which means 

fem(0) = 1 (12) 

if the meson wave function is properly normalized. 

In the calculation of the form factor at t ^ we start by using the 5^ functions to 
eliminate the integrals over the momenta q and k in the expression of T!p and over p and 
k in the expression of Z} 2 \ After performing the traces over 7 matrices we get 



T (l) = l-'M 



A* 



/ de p dp z d(p p d 4 Q — 5(e p + e q + Q - E) 5(e p - e k ) 
VI J Se p eke q 
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x <Pi(p, q\ Q) v>f{k, q; Q){q^ t + 2P ■ Q (k^ - 

+ (k^ + Pil )[{E-Q f + \t~Q 2 -{m l -m 2 f}} (13) 
and a similar result for T^ 2 \ Next, by writing 

~ Cp) = ^{jp z - P) (14) 



and 



1 

Ye 



S(e p + e q + Q - E) = 

-S I cos0 p - (E-Q ) 2 -2e p (E-Q o) + P 2 -Q 2 + m 2 - m 2 



2 Pt Qt \ v 2p T Q 7 



we perform the integrals over p z and (fi q in eq. (0). 

if' 



Then the term Tn 1 ^ becomes: 



(i) 



2ttk\ 1 
~VT~ V4M 2 - t 

{2e p 



2ptQt\JI - cos 2 
+(E-Q ) t}. 



d A Q 



{E-Q 



de p ipi(p,q; Q) ff(k,q; Q) 
, 1 
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t - Q - (mi - rn 2 ) 



(16) 



The integration limits over e p result from the kinematical constraints e 2 > m 2 + P 2 
and cos 2 p < 1 which give: 



where 



e pM 



(E - Q )S + Q T JS 2 - 4{(E - Q ) 2 - Q 2 ](m 2 + P 2 ) 



2[{E - Q ) 2 - Q 2 T ] 



Max[0, 



(E - Q )S - Q T JS 2 - 4[(E - Q ) 2 - Q 2 ] (m? + P 2 ) 



2[(E-Q ) 2 -Q 



(17) 
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S = (E- Q ) 2 + P 2 -Q 2 + ml- m 2 2 . 

The term Tq® can be proccessed in the same manner, giving a similar expression. 

Using the above results it is possible to calculate the electromagnetic form factors 
for any momentum transfer, by choosing an appropriate function ip. In principle, the 
calculation does not imply any other approximations, but it is hard to believe that the 
multiple integral entering the expression of the form factor can be performed exactly. 
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The numerical results quoted in this paper have been obtained in the approximation 
\Q\ << \Qo\ in the meson rest frame, in agreement with the assumption we made about 
the signification of Q^. 

We used the particular Lorentz invariant distribution function tp defined as 



v(p> Q',Q) = N ex P 
and performed the approximation 



[P-Qf-M 2 Q 2 

(3 2 M 2 





\(P-Q)] 


exp 


Ma 



(19) 



/ 7r \3/2 yr 

<p(p, q; QMk, q- Q) » (3 3 N 2 -5^ (Q) exp 



2P 2 Q 2 



M 2 f3 2 



exp 



2EQ 
Ma 



(20) 



expected to be valid for a small parameter f3. 

The approximation (^0|) allows one to do immediately the integration over Q and e p 
in 7q leaving only the integral over Q to be performed. The integration limits flI7D 
generated by the kinematical constraints become now: 



(E - Q ) 2 + P 2 + m 2 -m 2 
2(E-Q ) 

(E-Q ) 2 + P 2 -m 2 + m 2 



>m\ + P 2 



> m 2 , 



2(E-Q ) 

leading to a single condition for the integral over Q in T (1 ' , namely: 
Qo < Qom = Min[0, \/m 2 + P 2 - m 2 - \/ ' P 2 + m\\ 



(21) 
(22) 



(23) 



Performing the same operation in Tq 2 ^ and using the normalization condition (|5|) to 
eliminate the constant N, we finally get: 



fem(t) 



X < Ki 



7T M 2 



JPT E 2 



Qomi 



dQo exp 



\4M 2 (3 2 / 



exp 



2EQ 

M(3 



(E-Q c 



1 



2e p ( (E - Qo)' 2 - ~t - (mi - m 2 ) 2 



-«2 



QoM2 



dQ exp 



X 



4M 2 /3 2 
1 



exp 



2#q 

M/3 



t{E - Q c 



(E - Qo) 



2e q ( (E - Qo) - ^ - (mi - m 2 )" ) +t(£-Q 



(24) 
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where t = -4P\ e„ = {E - Q % E l'^"' and 

J = f_JQ« 1 - ' g ~ ^ exp (2^) v/[(A/ - Q„P - mf - m|] 2 - 4m>|. (25) 

It is easy to see that f^ v,v ' (t) = 0, while f^ ,K °(t) ~ (m s — m^) and in principle it 
does not vanish. 

The expression (|24|) is, of course, valid for t ^ 0, but the infinite value one gets in the 
limit t — > seems to contradict the normalization of the electric charge (^) which has 
been demonstrated previously. 

This is a disturbing question which deserves a careful examination. Looking back, we 
remark that the contradiction comes from the evaluation of some 5 functions: 

5 (3) (p + q + Q ~ P)S(e p + e q + Q - E{P)) 

x5 {3) (k + q + Q' - P')5(e k + e q + Q - E(P'))8^(Q - Q') (26) 
which have been written as 

6®(f+ q + Q - P)5(e p + e g + Q - E(P)) 

x6®@- k - 2P)5(e p - e k )5^\Q - Q')5(Q - Q' ) (27) 
at t t^O, while at t=0 they have been written as 
5 {3 \p + q + Q)8(e p + e q + Q -M) 

x6®(p- k)5(Q -Q' -M + M')5^\Q - Q')±- [ S M ~ M '^dX Q (28) 

and the integral has been replaced by T because it was assumed that the uncertainty in 
the meson mass is much smaller than T~ x . 

The problem comes from the fact that T is finite and hence it is illegal to put 8(Qo — Q' ) 
in the expression of the vacuum expectation value (|j). This means that instead of (|27|) 
we ought to write 

S(e p + e q + Q - E(P))5(e p -e k + Q -Q' - E(P) + E'(P')) 

X( J(3)(£ + Q _ p)5®(f- k - 2P)— [ e i(E(P)~E'(P')~e p+ e k )X odXo (2Q) 

and perform the calculation with T finite by also taking into account the indetermination 
of the meson mass. In the present calculation we do not follow this line because it is very 
cumbersome. Instead of this, we use the charge normalization condition (|T2| ) in order to 
fix the parameter T, which is mainly the same thing. 

To this end we notice that T is the overlapping time of the complex systems repre- 
senting the initial and final mesons. Then, as resulting from a careful analysis of the 
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relations (|27|) and 



one must write T = ^ where 



M 4 



- 1 



is the 
can now 



relative velocity of the two mesons. This solves the problem and the limit t 
be freely performed in eg . (|24l) . 

By using different values for the cut-off parameters a and (3 we found that the charge 
radii increase with a and decrease when the parameter (3 increases. We also found that 

2 

the shape of the electromagnetic form factor /(£) depends on the ratio p = For p > 1 
the shape is exponential, leading to large values for the charge radius, while for p < 1 it 
changes and the radius can be as small as wanted. 

The dependence of the shape on the ratio p is illustrated by the plots of the pion form 
factor in Fig.l. 




0.05 0.10 0.15 0.20 

-t (GeV 2 /c 2 ) 



Figure 1: Plots of the pion electromagnetic form factor for the following values of the 
parameters: a = 0.3 M w ; (a) /3=0.02 GeV; (b) ,3=0.04 GeV; (c) /3=0.06 GeV. 



For the comparison with the experimental data we shall retain however only the cases 
with p > 1, which fit the expected growth of the form factor at time-like t and are also 
in agreement with our initial assumption Q 2 > 0. 

The quark masses used in the calculations have been determined together with the 
cut-off parameter a from the fit of the decay constants of pseudoscalar mesons [H 
take: m u = 7MeV, m d =10MeV, m s = 400MeV and a = 0.3 M, which give F w : F K 
Fb '■ Fd s = 130. : 160 : 254 : 144 : 386., in agreement with the experimental data. We 
note that the values of the quark masses are rather close to the values suggested by chiral 
symmetry scheme |TT|. 



We 
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Moreover, by using the normalization condition ([12]) we get in the charged pion case 
To ~ 10~ 22 s, which is in agreement with the low values for the quark masses we used. 

Taking now p m l.(/3 = 0.04 GeV in the pion case), we find the following values for the 
charge radii: r 2 ± = 1.6 fm 2 , r 2 K± =l.Q fm 2 , r^o=-0.17 fm 2 , which are much larger than the 
measured ones: r„-±=0.44 fm 2 , r K ±=0.29 fm 2 , r^-o=-0.054 fm 2 [O]. We notice, however, 



„2 



the negative sign of r K o in agreement with the experimental result. This shows that the 
contribution of the heavy quark is dominant. 

The values we obtained for the charge radii suggest that the approximation (^) is 
inadequate. A simple way to improve it is to replace the symmetry scheme based on 
the full Lorentz group with the symmetry under the collinear group which is equivalent 
with the flux tube model with frozen transverse degrees ||13|| . The longitudinal and the 
temporal degrees of freedom are then the only active and the multiple integral in eq. ([TBI) 
reduces to a simple one. By using less drastic cuts of the internal momenta we expect to 
obtain a slower decrease of the form factors and a better agreement with the experimental 
data. The work on this line is in progress. 
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